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Using an effective field theory we describe the low energy bosonic excitations in a three dimensional 
ultra-cold mixture of spin-1 bosons and spin-1/2 fermions. We establish an interesting fermionic 
excitation induced generic damping of the usual undamped long wavelength bosonic collective Gold- 
stone modes. Two states with bosons forming either a ferromagnetic or polar superfluid are studied. 
The linear dispersion of the bosonic Bogoliubov excitations is preserved with a renormalized sound 
velocity. For the polar superfluid we find both gapless modes (density and spin) are damped, whereas 
in the ferromagnetic superfluid we find the density (spin) mode is (not) damped. We find that this 
holds for any mixture of bosons and fermions that are coupled through at least a density-density 
interaction. In addition, we predict the existence of the Kohn anomaly in the bosonic excitation 
spectrum of Bose-Fermi mixtures. We discuss the implications of our many-body interaction results 
for experiments on Bose-Fermi mixtures. 
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The interplay of bosons and fermions is ubiquitous 
throughout physics, ranging from the interaction of light 
(i.e. photons) and matter (i.e. electrons) to the behav¬ 
ior of a simple metal in the ionic lattice background. 
In solid state physics the interaction of electrons with 
slowly moving phonons provides the necessary attractive 
electron-electron interaction to form Cooper pairs lead- 
iirg to superconductivity [l|. While on the other hand 
the effect of the electron Fermi surface on bosons appears 
naturally in phonon excitations in the form of the Kohn 
anomaly Q. Another well-know example of the solid 
state manifestation of fermion-boson interaction is the 
polaron formation in ionic insulators. In the description 



(relevant for 
^) the fermions 


of itinerant quantum phase transitions_ 
various strongly correlated materials 
are always coupled to a bosonic collective mode reflect¬ 
ing the ordering of the underlying Fermi gas. A coupled 
fermion-boson interacting many-body system is thus a 
fundamental paradigm in condensed matter physics lead¬ 
ing to a large number of interesting phenomena. 

Coupled fermion-boson systems have also become of 
interest recently in ultra-cold atomic systems. In the 
context of cold atoms, a variety of Bose-Fermi mixtures, 
e.g., ®Li-7L^,|I^, 40K-87Rb [llHi3, ®Li-i33Cs 0,113, 
®Li-^^''’Yb [20j, have been prepared in experiments for dif¬ 
ferent purposes such as implementing symmthetic cool¬ 
ing |9l4ll|. studying molecule formation O, engineer¬ 
ing dipolar quantum simulators [l5| - exploring few-body 
physics IIThIQII or looking for interesting collective ex¬ 
citations [20|. In a recent experiment [20j of particular 
relevance to our theory to be presented in the current 
work, laser spectroscopy was used to study the effect of 
the fermions on the bosonic excitation spectra of ®Li- 
atomic mixtures. 

In the absence of fermions, the low energy excitations 
in Bose-Einstein condensates are well described by Bo¬ 
goliubov theory and it is well understood that the Bogoli¬ 


ubov quasiparticles (BQs) are damped through higher or¬ 


der interactions in the form of Beliaev |21l - l26l | and Lan¬ 
dau [23 damping. Due to destructive quantum inter¬ 
ference Landau and Beliaev process are suppressed [13 
at low momentum and low energies, thus making the 
long wavelength collective mode a well defined undamped 
bosonic excitation; in fact, the long-wavelength damping 
goes as vanishing rapidly as g —>■ 0. An important 
question of fundamental interest, which has also become 
relevant in view of recent experiments [Slii , is, how¬ 
ever, still open in spite of extensive theoretical activity, 
namely, how the fermion-boson interaction (specifically 
the existence of the Fermi surface) affects the bosonic 
excitation spectrum in a Bose-Fermi cold atom mixture. 
We address this important question in the current work 
using field theoretic techniques, finding a generic fermion- 
induced damping of the long wavelength bosonic collec¬ 
tive modes. 

We start with a microscopic Hamiltonian for a Bose- 
Fermi mixture, and derive a low energy effective field the¬ 
ory based on a controlled perturbative expansion, from 
which the low energy bosonic excitations and damping 
effects are obtained. A mixture of a spin-1 Bose gas and 
a spin-1/2 Fermi gas is considered and spin SU(2) sym¬ 
metry is assumed for both theoretical simplicity and rel¬ 
evance to experimental systems in the absence of mag¬ 
netic fields. We show quite generally that the linearly 
dispersive BQs of a bosonic superfluid interacting with 
fermions become damped due to Fermi surface effects, 
with a damping rate 


7(q)/;i = T>|q|, 


( 1 ) 


where V is dependent on the microscopic details of the 
system. Moreover, at large momentum (g Ri kp) we show 
that the Kohn anomaly appears as a kink in the BQ 
exception spectrum. The linear momentum dependence 
in Eq. m is drastically different from pure boson sys- 
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or Bose-Fermi mixtures with a Fermi 
2^, . On the other hand, 


terns 

surface that is gapped out 
the functional form of BQ linear dispersion at low energy 
Wsiql remains intact, with a renormalized sound velocity 
Vs- For spin wave excitations, we find damping in a polar 
(P) ground state whereas for a ferromagnetic (FM) state 
the spin waves are found to be undamped. 

We begin with the model fermion-boson interacting 
Hamiltonian, H = J d?r {'Hb +'Hf +'Hbf), with 


Hb = 


+ 


2mB 




T , . 

2 ^CL^a' 


J- cl)l (j) ,<J) 

^ 2 a 


TLf = 


2m F 


Vc^VCcr, 


( 2 ) 

( 3 ) 

( 4 ) 


TLbf = Ubf^I^c 


ctco 


J^lTab^^b 


(5) 


where repeated indices are summed over, T denotes 
a vector of spin-1 matrices, and cr denotes a vector 
of Pauli matrices. The operator dJa destroys a boson 
in the mz = a state and Ca destroys a fermion with 
spin a. The interactions Uo and U 2 are related to the 
scattering lengths of each hyperfine state of the bosons 
through Uo = {go + 2gf )/3 and U 2 = (gf - g^)/3, 
with g/ = ATrh^ajJrnB for a scattering length a/ in hy- 
perhne state / [3l|. Whereas the Bose-Fermi interac¬ 
tions Ubf and J are given hy Ubf = {Qifi + ‘^ 93 / 2 )/^ 
and J = 2(gf/f - gf/f)/3, with gfj^ = 2'nV?aF,,J'mBF 
as the s-wave scattering leng th with total spin Ftot and 
mBF the reduced mass [S^. We focus on a repulsive 
density-density Bose-Fermi interaction Ubf > 0 and fer¬ 
romagnetic spin-spin interaction J < 0. We assume the 
fermions to be non-interacting, which is valid for bare re¬ 
pulsive interactions because they are strongly irrelevant 
in the low energy limit [H, . 

In the absence of the Fermi gas, the spin-1 Bose gas 


can become [31j, 135| either a FM superfluid for U 2 < 0 


or a P superfluid for U 2 > 0, and we consider both situ¬ 
ations theoretically. In the FM phase, the ground state 
breaks both the U(l) and SU(2) symmetry of the Hamil¬ 
tonian and as a result hosts two distinct types of Gold- 
stone modes corresponding to gapless density excitations 
which are of the BQ form that go as |q| and of the fer¬ 
romagnetic spin wave form that go as q^. The P super¬ 
fluid also breaks the U(I) and SU(2) symmetry and also 
hosts two distinct sets of gapless density and spin wave 
excitations, however in this case they both take on the 
linear-in-g BQ form. 

We are interested in the low energy theory of the sys¬ 
tem and derive the corresponding effective action for the 
gas within a path integral framework (^ . The model 
Hamiltonian corresponds to an action in the grand canon¬ 
ical ensemble ^ = f dT(^H{T) + f d^r[‘^>l^T‘^^a + cl^TCcr — 

= Sb + Sf + Sbf- We have in- 




- 
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FIG. 1: Diagram (a) is referred to as the “tadpole” and (b) 
is the spin dependent particle-hole “bubble”. A solid line de¬ 
notes a fermion Green function with spin a, a solid rectangle 
denotes the Bose-Fermi interaction {Ubf or J), and the wavy 
line denotes the bosonic field cj)a in one of the niz = 1, 0, —1 
states. 


troduced the Bose gas action S'b, the Fermi gas action 
Sf, and their mutual interaction Sbf- To derive the 
low energy effective action of the bosonic superfluid we 
first expand the Bose operators about the mean field 
ground state -|-((/)i, (j>o, which together 

with the mean held value of gs gives rise to an effec¬ 
tive action for the (j) degrees of freedom (see 

the Supplemental Material). As the Fermi gas is non¬ 
interacting we integrate them out and determine the ef¬ 
fective action of the huctuations (/) (ignoring any con¬ 
stants) .?[(/)] = — Trlog(l — V{(j))Go)- The trace 

is over spin, space, and imaginary time indices and we 
have dehned the two by two matrices V {(p) and the Green 
function of the fermions Go, which depend on the super- 
huid state. As V{(p) is composed entirely of fluctuations 
of the Bose gas we can expand the logarithm to obtain 
the low energy effective action to quadratic order 

Ses[cp] = SbM^] + Ty{V{(P)Go) + ^Tr (^{V{cP)Gof) . 

( 6 ) 

The hrst term after SB,eS corresponds to the tadpole di¬ 
agram and the second to the spinful particle-hole bubble 
in Figs. [11(a) and (b) respectively. 

Ferromagnetic Ground State. Focusing on the ferro¬ 
magnetic case U 2 < 0 {{^\^pT^mf) = poz, where go is 
the boson density), the mean field expectation value of 
the Bose operators can be written as a three component 
spinor, ^J^f = (v^i 0, 0). The chemical potential of 
the bosons at the mean field level is = g^ Po + JMz + 

Ubfbf, where Mz = {{n^) - (n;))/2 and uf = 
with Her = cJ-Ccr. Expanding the Bose operators about 
their mean field ground state results in the fermions see¬ 
ing an effective magnetic field go J in the z direction and 
it breaks the up-down spin symmetry of the Fermi gas. 
This gives rise to the non-interacting Green function, 
Gq^^, = -S^^a'{dT-fi'^/{‘2mF)'^'^ - p-F + Jpocr/2), where 
CT = ±1 for t and and we have shifted the fermionic 
chemical potential pF = Pf — PqUbf- The matrix Vfm 
is given by 

Vfm{(P)= (^J{yJpQ/2(pl + ]^S'^)(j++\i.c^ (7) 

+ UBF{\fpo[(p\ + (pi] + n^)l -t- —{\fpo[(p\ + (pi] + 

where 1 denotes the identity matrix, Ui are the Pauli 
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matrices, (7± = ± iay)l2^ and we have introduced 

= Ea and 

Evaluating the trace for the particle-hole bubble leads 
to XapiQ) = f dkGoa{k + q)Gof 3 {k), we are using the 
shorthand notation q = {wn, q) [k = {iujn, k)] for bosonic 
[fermionic] Matsubara frequency, and the integral f dq = 
f "t’tie Green function in momentum space 

is Gocr(k, — Ckcr + p-F, with a spin de¬ 

pendent dispersion Ckcr = fi^k^/(2mF) + Jpo<7/2. For 
the spin diagonal case a = B this reduces to the well 
known Lindhard function [37|, which at sufficiently low 
temperature in the low energy limit, \q\/kF^ <C 1 and 
i/„/(ftn_Fcr|q|) <C 1 (where vpa is the Fermi velocity of 
the spin a electrons), becomes 


= -Pa 


TT Wn\ 

2 hvFcr\(l\ 


1 

3 



( 8 ) 

and we have defined the density of states for spin cr, 
? 7 cr = mpkpa/{2-K'^fi^). We stress the approximation 
i/„/(ftu_Fcr|q|) <C 1 is consistent while considering Bose- 
Fermi mixtures with z;s q “C vp^ where Vgfi is the sound 
velocity of the BQ excitations in the absence of the Fermi 
gas. Physically, this ensures that the BQs excite a suf¬ 
ficient amount of particle-hole pairs, by passing through 
their excitation continuum as shown in Fig. (d). For 
the case a ^ a', in the low frequency, low momen¬ 
tum limit we can always treat h‘^\q\'^/{2mp) <C |J|po, 
hvpa\ci\ ^ \ J\pq, and |t'„| <C | J|po. Within this approxi¬ 
mation we find 




q,Wn) = -7 ((«<t) - {n, 

PoJ 


,, , _ IVn 


'h?np 
- 2mp 


+ i[hvpgf-{ns) - {hvpgf{n^)) 
^PqJ 


(9) 




With all of these results in hand we can now de¬ 
termine the effective action to quadratic order. Writ¬ 
ing the action in separate parts we have = 

/ dq -I- I where 


>^ 4>1 = '/’!('?) {-i^n + boq^) (j3i{q) 


( 10 ) 


(A™ + \ct>\{q) + M-q) 


+-B^^|qp)(/)J(g) 0 o(( 7 ), and = 

[M - iFn + bo\q\'^) (t)Liiq)4>-i{q), where bo = h?/2mB. 
In order to simplify the presentation we have introduced 
the constants AfSfand B™, whose 
functional form is not particularly relevant for the present 
discussion and are explicitly given in the Supplemen¬ 
tal Material. Now that we have the Green function for 
each spin state from the effective action, we can deter¬ 
mine the excitation of each mode from the poles [^ . 
A few remarks are in order: The (j>-i mode is gapped 


with a value A4 = 2{po\U2\ + M^jJl). In contrast, the 
spin wave excitations that correspond to (po, are given by 
ujfjf = /{2ml) where the coefficient is altered from 

/{2m b), and the spin excitations acquire a renormal¬ 
ized effective mass fP'/{2ml) — Bq^/Aq^, which is not 
a function of Ubf- Therefore, the renormalization of M 
and m* is due entirely to the paramagnon excitations 
of the Fermi gas through the coupling J. The density 
modes {cpi +(j)\) acquire damping through the additional 
contribution of li^nl/lql which arises due to the particle- 
hole excitations of the Fermi gas (see Eq. (ITT)) below). It 
is useful to note that the constants Af ^ and B™ are 
both quadratic functions of Ubp and J. 

Polar Ground State: We now discuss the polar ground 
state of the bosons when U 2 > 0, ((<I)]^^pT$MF) = 0), 
with a mean field expectation ^]/[p = (0, ^/pQ, 0). Now, 
the up down spin symmetry of the Fermi gas remains 
intact. The chemical potential for the bosons takes 
the form p// = UqPq + UBpnp. The non-interacting 
Green function is now spin independent, which becomes 
= -da,a'{dr - / (2mB)V^ - pp). The matrix Vp 

is defined as 

Vp{(j)) = Ubf{^/po['Po 'Po] + ncj,)l + JS^az/2 
+ Po/‘^[P-i + 4>i] + i^S^ )(T+ -I- h.c.^ . (11) 

At this point it is useful to compare this with the FM 
case. From the effective bosonic action SB,es[<P] in the P 
case we know that the density mode is related to Sn{q)'^ = 
{Po{q),po{—q)) while the spin excitations are given by 
SS{q)^ = {p\{q),p-i{—q)). As a result, it is quite natural 
that the density excitations only couple through Ubf and 
the spin excitations through J. This is quite different 
from the case of the FM superfluid where the BQs are a 
combination of density and spin along the z direction. 

Following the same steps as before, we compute the 
trace in Eq. only now the equal spin particle-hole 
bubble in Eq. ([8|) comes into play Due to the spin sym¬ 
metry the Fermi gas parameters are now a independent. 
The effective action to quadratic order for the polar case 
is given by, 5^[(?(>] = J dq {^Csn + Css) , where Csn = 
dn{q)''^Gsn{q)~^Sn{q) and Css = 6S{q))Gss{q)~^5S{q). 
The bosonic Green functions can be written in terms 
of the free bosonic Green function and the self energy 
using Dyson’s equation Ga(*r'„,q)“^ = G°(iF„,q)“^ — 
Ca{ivm q), which are two by two matrices G/^{iVn, q)~^ = 
—ii'n<7z -I- boq^l and 


Ca{iFu,q) = - (^Af -I--f Sf |qp^ (l-kcr„)(12) 

We give the explicit forms of A//, A//, and B// in the 
Supplemental Material. In this case, the 6 n constants 
are quadratic functions oi Ubf and do not depend on 
J while the 5S constants are quadratic functions of J 
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FIG. 2: Energy spectrum and damping rate of collective 
modes from solving Eq. m for various bosonic densities po- 
In this plot the exact Lindhard function [s^] is included fo¬ 
cusing on the density mode (Sn) for the P case with ^^Na for 
the bosons, the Fermi gas it taken to be ®Li with a density 
of n_F = 10^^cm“®, and focusing on Ubf = Uo, and J = U 2 - 
(a) The dispersion Esn{<l), (b) the damping rate 75 n(q), (c) 
the difference in dispersions with and without the Fermi gas 
5Esn{<l) = Esn,o{<l) — Esn{<l), displaying the Kohn anomaly 
near q/kp ~ 1.5, and (d) the Bose gas dispersion in the ab¬ 
sence of the Fermi gas i?^„,o(q) = + b^q^, all as a 

function of q/kp- In (d) the black solid lines bound the re¬ 
gion of the particle-hole excitation continuum and the black 
dashed line marks where the imaginary part of the Lindhard 
function changes its form. 


and are independent of Ubf- Interestingly, in the polar 
case we find both BQ modes to be damped through the 
appearance of li/nl/lql. We remark that the density mode 
for the FM case in Eq. (ITUl) can also be converted to the 
above form. 

BQ Excitation Spectrum and Damping: We have de¬ 
rived the effective field theory for both the ferromagnetic 
and polar bosonic ground states interacting with a spin- 
1/2 Fermi gas. To diagonalize the effective action using 
the Bogoliubov transformation we find rtk and Uk have 
to be functions of Pn- To determine the dispersion and 
damping rate we find the poles of the bosonic Green func¬ 
tion Ga (with a = (l>i,6n,SS) [s^, after analytically 

continuing the Matsubara frequency to real frequency 
ivn -P- Lo fO’*', we solve the algebraic equation 

det Ga(w, q)"^ = 0, (13) 

to find the poles at fiw = ^(q) —i 7 (q), where 7 (q)//i > 0 
is the damping rate. We obtain two of our main results 

7 a(q) = boAa\q\, Ea{q) = + &2q4, (14) 

(see Fig. E]) and we find the Bogoliubov form remains with 
the renormalized parameters uf = bo{2Aa — bgA^), and 
ba = bq {bo + 2Ba) ■ The low energy theory can also be 
extended to finite temperature, which yields 7 a(q, T) oc 


[1 -E tanh{flF/kBT)]\q\ and the linear momentum depen¬ 
dence survives. Solving Eq. (fT3l) using the exact Lind¬ 
hard function gives results indistinguishable from the 
analytical results as shown in Fig.[2](b) up to g ~ 0{kp). 
Thus, all BQ modes present in a spin-1 Bose gases cou¬ 
pled to a Fermi gas become damped at low energies and 
are thus no longer true collective modes even at long 
wavelengths! 

Several comments about our main results shown in Eq. 
(fT4)l and Fig. [2] are in order: (i) Our results are valid to 
all orders in the Bose-Fermi coupling provided Vsp <C vp 
since all higher-order vertex corrections in the self-energy 
are negligible by virtue of Migdal’s theorem 3^. (ii) 
The fermion-induced bosonic damping being linear in q 
strongly dominates any intrinsic bosonic Beliaev damp¬ 
ing at long wavelength [ 2 ^. (hi) Eq. da implies that 
the bosonic collective mode frequency vanishes without 
becoming overdamped if the damping becomes compara¬ 
ble or larger than the mode energy itself since both the 
frequency and the damping go linear in wave number. 

(iv) At large wave numbers, q > kp, 7a (q) —>■ 0 because 
the imaginary part of the Lindhard function [^1 itself 
vanishes by energy-momentum conservation and the BQs 
can no longer excite particle-hole pairs [see Fig. [2] (d)] 

(v) The spinor structure of the problem is not essential 
to find this Fermi gas induced damping, therefore our re¬ 
sults are valid for any condensed spin-S' Bose gas coupled 
to a Fermi gas (either spinful or spinless) through at least 
a density-density interaction UBpnBnp. 

The Fermi surface has effectively been imprinted upon 
the excitation spectrum of the Bose gas in the form of 
the Kohn anomaly. As this is a large momentum effect, 
it is not captured by our low energy theory but comes 
directly out of the numerical solution using the exact 
Lindhard function. However, this effect is weak for the 
atomic mixtures we are considering and is not visible in 
ifa(q) itself, only the difference SEa{q) = Eo{q) — Ea{q), 
displays the anomaly as seen in Figs. [2] (a) and (c). 

Our theoretical predictions can be tested in ultra cold 
Bose-Fermi mixtures with weakly interac ting fermions, 
through two-photon Bragg spectroscopy [28( that cou¬ 
ples to the bosonic species. We expect the damping will 
give rise to a broadened Bragg line, i.e. an intrinsic line 
width, in the presence of fermions that will depend ex¬ 
plicitly on energy-momentum as shown in Eq. (1141) . Due 
to the bosonic damping persisting to sufficiently low mo¬ 
mentum, we also expect that the presence of the Fermi 
gas can dephase the bosonic superfluid, which has possi¬ 
bly already been observed in Ref. 

In the presence of a shallow harmonic trap as in most 
atomic experiments, we expect our results to be largely 
unaffected for large clouds of atoms. Even though the 
finite size of the trap gives an infrared cut off to the 
theory, this energy scale is sufficiently small such that 
linear-in-g BQ dispersion can still be observed in exper¬ 
iments [39|. In this limit, the physics is well described 
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by the Thomas-Fermi and local density approximations, 
where the BQs will still be linear-in-g at low momentum 
(with a renormalized sound velocity) [si, H, 4^ and the 
Fermi gas will acquire a spatially dependent Fermi ve¬ 
locity. In spite of this, close to the center of the trap, 
the form of particle-hole excitations will remain the same 
with a renormalized vp- However, we would like to men¬ 
tion that these inhomogenous effects can in principle be 
removed by considering a uniform potential trap, which 
has recently been achieved, where the spatial inhomo¬ 
geneity is minimized 41|. 


We have studied a spin-1 Bose gas coupled to a spin- 
1/2 Fermi gas in three dimensions. We have shown 
BQ excitations are damped at long wavelength as a re¬ 
sult of particle-hole excitations of the Fermi gas while 
the functional form of the BQ dispersion is unchanged. 
We have argued this phenomenon should apply to Bose- 
Fermi mixtures in general independent of the spin struc¬ 
ture of either species. 


Acknowledgements We would like to thank Stefan Natu 
and Brian J. DeSalvo for useful discussions. This work 
is supported by JQI-NSF-PFC and ARO-MURI. 


[1] J. Bardeen, L. Cooper, and J. Schrieffer, Phys. Rev. 108, 
1175 (1957). 

[2] W. Kohn, Phys. Rev. Lett. 2, 393 (1959). 

[3] J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 

[4] A. J. Millis, Phys. Rev. B 48, 7183 (1993). 

[5] T. Ishiguro and K. Yamaji, Organic superconductors 
(New York, NY (USA); Springer-Verlag New York Inc., 
1990). 

[6] A. C. Hewson, The Kondo problem to heavy fermions, 2 
(Cambridge university press, 1997). 

[7] N. Armitage, P. Fournier, and R. Greene, Rev. Mod. 
Phys. 82, 2421 (2010). 

[8] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, 
Journal of the American Chemical Society 130, 3296 
(2008). 

[9] A. G. Truscott, K. E. Strecker, W. I. McAlexander, G. B. 
Partridge, and R. G. Hulet, Science 291, 2570 (2001). 

[10] F. Schreck, L. Khaykovich, K. L. Corwin, G. Ferrari, 
T. Bourdel, J. Cubizolles, and C. Salomon, Phys. Rev. 
Lett. 87, 080403 (2001). 

[11] G. Roati, F. Riboli, G. Modugno, and M. Inguscio, Phys. 
Rev. Lett. 89, 150403 (2002). 

[12] J. J. Zirbel, K.-K. Ni, S. Ospelkaus, J. P. D’Incao, C. E. 
Wieman, J. Ye, and D. S. Jin, Phys. Rev. Lett. 100, 
143201 (2008). 

[13] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe’er, 
B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, 
D. S. Jin, and J. Ye, Science 322, 231 (2008). 

[14] S. Ospelkaus, K.-K. Ni, D. Wang, M. H. G. de Miranda, 
B. Neyenhuis, G. Qumner, P. S. Julienne, J. L. Bohn, 
D. S. Jin, and J. Ye, Science 327, 853 (2010). 

[15] A. Chotia, B. Neyenhuis, S. A. Moses, B. Yan, J. P. 
Covey, M. Foss-Feig, A. M. Rey, D. S. Jin, and J. Ye, 
Phys. Rev. Lett. 108, 080405 (2012). 


[16] K. Gunter, T. Stoferle, H. Moritz, M. Kohl, and 
T. Esslinger, Phys. Rev. Lett. 96, 180402 (2006). 

[17] R. S. Bloom, M.-G. Hu, T. D. Cumby, and D. S. Jin, 
Phys. Rev. Lett. Ill, 105301 (2013). 

[18] S.-K. Tung, G. Parker, J. Johansen, G. Ghin, Y. Wang, 
and P. S. Julienne, Phys. Rev. A 87, 010702 (2013). 

[19] S.-K. Tung, K. Jimenez-Garcia, J. Johansen, G. V. 
Parker, and G. Chin, Phys. Rev. Lett. 113, 240402 
(2014). 

[20] H. Hara, H. Konishi, S. Nakajima, Y. Takasu, and 
Y. Takahashi, Journal of the Physical Society of Japan 
83, 014003 (2014). 

[21] S. T. Beliaev, Sov. Phys. JETP 7, 289 (1958). 

[22] P. C. Hohenberg and P. C. Martin, Ann. Phys. (N.Y.) 
34, 291 (1965). 

[23] V. N. Popov, Theor. Math. Phys. 11, 478 (1972). 

[24] W. V. Liu, Physical Review Letters 79, 4056 (1997). 

[25] S. Giorgini, Phys. Rev. A 57, 2949 (1998). 

[26] S. S. Natu and S. Das Sarma, Phys. Rev. A 88, 031604 
(2013). 

[27] L. Pitaevskii and S. Stringari, Physics Letters A 235, 398 
(1997). 

[28] R. Ozeri, N. Katz, J. Steinhauer, and N. Davidson, Rev. 
Mod. Phys. 77, 187 (2005). 

[29] 1. Ferrier-Barbut, M. Delehaye, S. Laurent, A. T. Grier, 
M. Pierce, B. S. Rem, F. Chevy, and C. Salomon, Science 
345, 1035 (2014). 

[30] W. Zheng and H. Zhai, ArXiv e-prints (2014), 1408.6419. 

[31] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

[32] Z.-F. Xu, X. Li, P. Zoller, and W. V. Liu, arXiv preprint 
arXiv: 1407.3720 (2014). 

[33] R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 

[34] J. Polchinski, ArXiv High Energy Physics - Theory e- 
prints (1992), hep-th/9210046. 

[35] T. Ohmi and K. Machida, Journal of the Physical Society 
of Japan 67, 1822 (1998). 

[36] J. W. Negele and H. Orland, Quantum many-particle sys¬ 
tems, vol. 200 (Addison-Wesley New York, 1988). 

[37] J. Lindhard, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 
28, 8 (1954). 

[38] A. B. Migdal, Sov. Phys. JETP 7, 996 (1958). 

[39] J. Steinhauer, R. Ozeri, N. Katz, and N. Davidson, Phys. 
Rev. Lett. 88, 120407 (2002). 

[40] F. Zambelli, L. Pitaevskii, D. M. Stamper-Kurn, and 
S. Stringari, Phys. Rev. A 61, 063608 (2000). 

[41] A. L. Gaunt, T. F. Schmidutz, 1. Gotlibovych, R. P. 
Smith, and Z. Hadzibabic, Phys. Rev. Lett. 110, 200406 
(2013). 






6 


Supplemental Material 

In the supplemental material we give the explicit expressions for Sb^bS and the equations defining the constants 
that we have introduced in the main text in terms of the parameters of the model. For the ferromagnetic case, the 
effective bosonic action including the bosonic chemical potential is 


— J drcPr f {JMz + UBFnp) 
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whereas for the polar case, the effective bosonic action takes the form 

‘S's,eff[^] = I + “^(^0 +<Ao)(4 +«^o) + ^2^0(4 + ?^'l) I 
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To simplify the presentation of the main text we have defined the following constants for the ferromagnetic case 
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For the polar case we have introduced the constants 
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